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| Abstract. A general prescription for the treatment of constrained quantum motion 

is outlined. We consider in particular constraints defined by algebraic submanifolds 
of the quantum state space. The resulting formalism is applied to obtain solutions to 
' j—i ' the constrained dynamics of systems of multiple spin-i particles. When the motion 

O ■ is constrained to a certain product space containing all of the energy eigenstates, 

, 1 * , the dynamics thus obtained are quasi-unitary in the sense that the equations of 

motion take a form identical to that of unitary motion, but with different boundary 
conditions. When the constrained subspace is a product space of disentangled states, 
Q" 1 the associated motion is more intricate. Nevertheless, the equations of motion satisfied 

by the dynamical variables are obtained in closed form. 
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Practical implementations of quantum algorithms in quantum information 
technology can be hampered by constraints. The formulation of a tractable theory of the 
properties of constrained unitary motions is therefore of interest and importance. This 
paper provides a framework for dealing with certain classes of constrained quantum 
motion, for which differential equations governing the constrained dynamics can be 
obtained explicitly. The results are applied to constrain unitary motions to algebraic 
subspaces of quantum state spaces. 

The structure and the findings of the paper can be summarised as follows. We begin 
by reviewing the Hamiltonian formulation of standard quantum mechanics. This highly 
effective way of looking at quantum mechanics has been investigated by a number of 
authors (see [IllSlEllIlElEllTJIEllH] and references cited therein) . We demonstrate the 
existence of a particularly simple choice of canonically conjugate coordinates {Pi}) 
for the space of pure quantum states with the property that the Schrodinger equation 
can be written in Hamiltonian form: 

a dfffap) dH(q,p) 

qi ~ d Pl V% ~ dq t ■ (1) 

The coordinates ({<&}, {Pi}) appearing here are defined by the expansion coefficients of 
the normalised state vector of a non-degenerate n-level system in terms of the energy 
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eigenstates {|-Ej)}i=i,2,...,n according to the scheme 

/ n-l N 

y /p i e-^\E i )+ (l-^Pi) \E n ), (2) 



n— 1 / n-l 

i=i \ t=i 



where the function 



is given by the expectation of the Hamiltonian operator in the state defined by (J2j). We 
choose the overall phase such that the coefficient of the nth energy eigenstate \E n ) is 
real. Since in the energy basis we have 

n 

H = ^E i \E i )(E i \, (4) 
i=i 

it follows that the Hamiltonian function is given by 

n-l 

H(q, p) = E n + uJiPi, (5) 

i=l 

where u)i = E{ — E n . Note that H(q,p) is independent of {qi} and is linear in {pi}. 
Substitution of (jSJ) in (pP) shows that the solution to the Schrodinger equation 

i^ t ) = m) (6) 

is given in terms of the canonical coordinates by 

Qi(t) = ft(0) +wib and Pi(t)=Pi(0), (7) 

which should also be evident from (T5]). For simplicity we shall consider Hamiltonians 
having nondegenerate eigenvalues, although the formalism can be applied to degenerate 
systems by use of the Liiders projection postulate [TO] . Specifically, the state \ip) is still 
expressible in the form (j2j), but we make the replacement n —>■ d, where d < n is the 
number of distinct energy eigenvalues, and \Ej) = U.^), where flj is the projection 
operator onto the eigenspace associated with the eigenvalue E iy given by 

di 

f[ i = Y,\Ei,j){E i ,j\. (8) 

3=1 

Here dj is the dimension of the Hilbert subspace associated with the eigenvalue E{, 
and \Ei,j) (j = l,...,di) constitute an orthonormal basis for that subspace. The 
Hamiltonian for a general system is hence given by 
d 

H = J2 E A, (9) 

and one sees that the form of (jSJ) remains unchanged, except that n is replaced by d, 
and H(q,p) is independent of the remaining phase space degrees of freedom. 

The objective of this paper is to introduce a framework for treating certain classes 
of constrained unitary motion. Our approach is aligned closely with that of Dirac's 
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theory of constraints in classical mechanics [III [12] • The idea that Dirac's methodology 
might be applied to investigate constrained quantum motion was proposed recently 
by Buric [13] to determine the dynamics of a pair of spin-| particles constrained to a 
special surface of product states containing all of the energy eigenstates. An alternative 
approach to quantum constraints is considered in [14J. 

An elementary way of enforcing the constraints is to introduce Lagrange multipliers. 
In some circumstances the Lagrange multipliers can be determined explicitly and 
eliminated from the equations of motion. In this paper we shall be considering such 
cases. Several examples are investigated, including one for which the motion of a pair 
of spin- 1 particles is constrained to the hypersurface of disentangled states. 

Constrained unitary motions are in general nonunitary, and correspond to nonlinear 
evolutions. However, unlike the general nonlinear dynamics of the Mielnik-Kibble- 
Weinberg (MKW) framework j3[ HI [15], the nonlinearities resulting from the class of 
constraints considered in the present investigation are of a milder form. That is, while 
in the MKW theory one considers a general Hamiltonian H(q,p) that is distinct from 
([5]), in the present context the 'linear' Hamiltonian (Q remains unchanged, and the 
nonlinearity arises from a modification of the symplectic structure, or equivalently on 
account of the nonlinearity of the constraint surface. 

In the first example we consider a system of two spin-| particles. The constraint 
surface is defined by the product space of a pair of two-dimensional Hilbert spaces, 
and we assume that the product space contains the energy eigenstates. An initial 
state that lies on the constraint surface is thus obliged to remain on this surface. We 
shall obtain the trajectories of the unitary evolutions subject to this constraint. The 
analysis can be extended and applied to n spin-| particles constrained to a special 
product space containing the energy eigenstates. We find that this constraint is satisfied 
under unitary evolution without constraint if the spectrum of the trace-free part of the 
Hamiltonian takes the form {Ei} = {e±, e 2 , • • • , e n / 2 , — e„/ 2 , — e„/ 2 -i, • • • , — ei}, where n 
is the number of eigenstates of the system. For a generic Hamiltonian we derive and 
solve the constrained equations of motion explicitly in the case of two and three spin-| 
particle systems. Surprisingly, the resulting dynamics turn out to be quasi-unitary in 
the sense that the amplitudes {p^ remain constant while the relative phases {q{\ evolve 
linearly. This property appears to be generic even for systems with more particles. We 
then consider the motion of two spin-| particles constrained to remain on the quadric 
corresponding to the subspace of disentangled states. 

Let us now begin with a summary of the Hamiltonian formulation of quantum 
mechanics. We consider a complex Hilbert space 7i n of dimension n, a typical element 
of which is denoted If F represents an observable, its expectation with respect to 
l^) is F = (ip\F\ip)/(i/)\ip), which is invariant under the transformation — > 
A G C — {0}. The vector thus carries a redundant complex degree of freedom. 
We therefore construct the space of rays through the origin of 7i n by the identification 
l^) ~ A|t/>). The result is the projective Hilbert space V n ~ x of dimension n — 1. We 
can view V n ~ x as a real even-dimensional manifold T, and let {x a } a =i,2,...,2ri-2 denote a 
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typical point of T. One distinguishing feature of T is that it is equipped with a symplectic 
structure Q ab (= —Q ba ) such that the Schrodinger equation ([6]) can be expressed in the 
Hamiltonian form 

± a = Q ab V b H. (10) 

Here the function H(x) = (ifj(x)\H\ip(x)) / (ip(x)\tlj(x)) denotes the expectation of the 
operator H in the pure state \ip{x)) corresponding to the point x G V . One can regard 
T as a bona fide quantum analogue of the phase space in classical mechanics. 

Now suppose that we have a family of constraints on the motion of the system in 
T expressed in the form 

$ a (a;) = 0, (11) 

where a = 1, 2, . . . , N. There are two different types of constraints that arise naturally 
in the quantum context, corresponding to what one might call 'algebraic' and 'real' 
constraints. In the algebraic case the motion is confined to an algebraic submanifold 
(or possibly a complex algebraic subvariety) of the original quantum state space V n ~ x . 
As a consequence, the constraint submanifold is of even real dimension — for which it 
follows that the number of constraints N is even in this case. This paper is primarily 
concerned with the algebraic case. Typical examples would include the situations where 
the constraint manifold 9Jt was an algebraic curve in V 2 (such as a conic or a elliptic 
cubic curve), or an algebraic curve in V 3 (such twisted cubic curve or an elliptic 
quartic curve), or an algebraic surface in V 3 (such as a quadric surface, or a cubic 
surface). The equations given by (TTTT) then define 9Jt locally. 

The other situation that is natural to consider in quantum theory is the case where 
the constraints are of the form 

where {F a } a =i,...,N denotes a collection of observables, {f a } a =i,...,N is a set of real 
numbers, and N need not be even. We shall consider the 'real' case elsewhere. 

In the algebraic case, the constraints can be enforced by the introduction of 
Lagrange multipliers {X a } a =i,2,...,N- Using the usual summation convention, the 
constrained equations of motion are 

x a = Q ab V b H + \ a Q ab V b <S> a . (13) 

To determine the Lagrange multipliers we analyse the relation $ a (x) = 0. From the 
chain rule we have $ Q = i; a V a ( I )Q! = 0. Substituting ( JT3l) in here, we find 

n ab V a $ a V b H + XpQ ab V a ^ a V b ^ = 0. (14) 

To solve (fT4"j) for X a let us define 

= n a6 V a $ a V fc $^. (15) 

In the case of real constraints, for which corresponds to a family of observables, 

u afS is the commutator of the observables $ a and We note that since £l ab = —tt ba 
we have co a/3 = —ui^ a . 
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If the matrix u al3 is nonsingular, then we can invert it. In that case, writing uo a p 
for the inverse of u al3 so uj^uj^ = 5 a ^, we can solve (HM for {X a } to obtain 

X a = u Pa n ab V a &V b H. (16) 
Substituting this in the right side of f|T3|) yields 

x a = Q ab V b H + wi 3a n ed V c $ f} V d HSV*V b $ a . (17) 
This can be simplified further by rearrangement of indices, after which we deduce that 

x a = Q ab V b H, (18) 
where & ab = Vl ab + A ab and 

A ab = Q ac Q bd uJ lS V c ^V d ^ s . (19) 

An important point to note is that A ab is by construction antisymmetric. Therefore, Q ab 
defines a modified symplectic structure. The constrained equation of motion (1181) thus 
takes on a form identical to ffTU]) . with the same Hamiltonian, but with the modified 
symplectic structure. 

The modified symplectic structure can be interpreted as playing the role of an 
induced symplectic structure on the constraint surface $ = 0. To see this we transvect 
Q ab with the vector V a $ a normal to the constraint surface to obtain 

Q ab V b $ a = Q ab V b $ a + Q^Q^sVc&'Vd&Vh^. (20) 
Using the antisymmetry of fl ab and the definition (fl5j) we find 

tt bd V d $ 5 V b $ a = -iu Sa . (21) 
Hence from uj 1 su 5a = 5" we deduce that 

VL ab V b $ a = (22) 

for all a. Therefore, Q ab annihilates all vectors normal to the constraint surface, and 
hence induces a symplectic structure on the constraint surface. 

Our procedure for dealing with a constrained unitary motion can be summarised as 
follows: (i) find a suitable choice of In — 2 real coordinates for representing the generic 
pure state (ii) calculate the symplectic structure fl ab in that coordinate system so 
that the unitary evolution is represented in the Hamiltonian form (|T0|) ; (iii) express 
the constraints ( JTTI) in terms of the given choice of coordinates; (iv) assuming that the 
constraints are such that the matrix uj af3 of ( JT5l) is invertible, calculate A ab according to 
f|T9|) and substitute the result into f|T8l) . In this way, dynamical equations for constrained 
unitary motion can be obtained, and one is left with the problem of solving a system of 
coupled differential equations. 

As we have indicated, there is a particular choice of coordinates on the space of 
pure states for which the analysis can be simplified in the form defined in (T5]), which 
might appropriately be called an 'action-angle' parametrisation (cf. [16]). In terms of 
these coordinates Q ab is given by 

n- = I ° I I - (23) 
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where O and 1 denote the (n — 1) x (n — 1) null matrix and identity matrix. As a 
consequence, the dynamical equations ffTUj) take the form ([I]). 

It is worth noting that while in classical mechanics phase space coordinates 
correspond to observables, in quantum mechanics only half of the phase space 
coordinates correspond to observables. Specifically, if we write IT = \Ei)(Ei\ for the 
observable corresponding to the projection operator onto the ith normalised energy 
eigenstate, then = (ip\fli\ijj) . Therefore the coordinates {pi}i=i t 2,...,n-i constitute a 
commuting family of observables. The conjugate variables {(?j}t=i,2,...,n-i, correspond to 
the relative phases, and do not represent observables in the conventional sense. 

Example 1. We now apply the formalism to specific examples. The first example 
is a system consisting of a pair of spin-| particles. For a generic Hamiltonian, we shall 
impose the constraint that under the dynamics the initial state of the system remains 
on a quadratic surface Q = V 1 x V 1 C T that contains the energy eigenstates. Such a 
constraint implies that the quantum state can be represented as a product state with 
respect to some choice of basis elements. The Hilbert space is four dimensional and a 
generic state can be expressed in the form 

|V>> = VPi^ iqi \Ei) + VP2e- iq2 \E 2 ) + VPse- iq3 \E 3 ) + (1 - Pl - p 2 - p 3 ) 1/2 \E 4 ). (24) 

The space of pure states is the three-dimensional space V 3 , in which sits the two- 
dimensional product space Q. The constraint for the state to remain on Q is therefore 
of an algebraic type — that is, \ip) must lie on the algebraic subspace Q. If we write 
f0t}i=i,...,4 for the coordinates of the Hilbert space vector \ip) = (ipi,ip 2 ,ip3,?p4), then a 
necessary and sufficient condition for to lie on Q is ^1^4 = ^2^3 [2]. Expressing the 
real and the imaginary parts of this condition in terms of the coordinates chosen in ( |24l ) 
we find that the constraint equations are given by 

y/pipi cos qi - jp^pi cos(g 2 + qs) = 
v/PiP4 sin qi - y/p 2 P3 sin(g 2 + q 3 ) = 0, 

where for brevity we have written p 4 = 1 — pi — p 2 — p%. If we divide the first equation in 
(|25|) by ^/p 2 pz cosgi and the second by ^/p 2 pz singx, and compare the results, we find 
that the constraint equations can be made separable: 

$ 1 = q\ - <?2 - qs , 0fi x 

x>2 /i \ ( 26 ) 

= Pl{l ~ Pl ~ p 2 ~ P3) - P2P3- 

Before we proceed to derive the constrained dynamical equations, we address the 
following question: What is the condition on H that will ensure that under unitary 
evolution an initial state that lies on Q remains on Q? The answer is obtained by 
substituting the solution ([7j) of the unitary motion into the constraints (f2"ET) . We thus 
obtain the condition that lo\ = uj 2 + Translated into the eigenvalues of H this 
condition is E x — E 2 = E 3 — E±. It follows that the trace- free part of the Hamiltonian 
must have the eigenvalue structure {ei, e 2 , — e 2 , —ex}. 

We now turn to the general case for which E\ — E 2 ^ E% — E±. As a consequence 
of the separable decomposition (1261) . the matrix uj al3 and its inverse uj a p are remarkably 
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Figure 1. A 'snapshot' of the vector field generated by (di,<fii) in Example 1. The 
parameters are chosen to be Ei = |, £? 2 = 1, E 3 = —2, E A = i, and 9 2 = ^n, thus 
giving 9i — and 4>i = \ — sin 2 \Q\. In this example, 0\ > in the northern 
hemisphere, 9\ = along the equator, and 6\ < in the southern hemisphere, 
resembling the nonlinear motion considered by Mielnik [17) . 



simple in this example: 
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(28) 



(29) 



where, writing p 4 = 1 — pi — p 2 — p% as before, we have 
' Pi - P4 Pi- Pi Pa- Pi ^ 

A = P1+P3 -Pi - P3 ~Pl ~ P3 

\Pl+P2 ~Pl ~ P2 ~Pl ~ P2 J 

Substituting these results into (fT8l) we obtain the following equations of motion: 

qx = (ujx -u 2 - ^3)(2pi +P2+ Ps) + (u>2 + u z ) 
q 2 = (ujt -u 2 - w 3 )(pi + p z ) + u 2 
Q3 = (wi -uj 2 - uj 3 ){Pi + Pz) + ^3 
Pi = 
P2 = 
P3 = 0. 

It should be evident that if the condition uj\ = oj 2 -\- u 3 holds, then (130!) reduces to the 
unitary case. It is interesting to observe that in spite of the fact that the evolution is 
no longer unitary we nevertheless have pi = and hence = constant for i — 1, 2, 3. In 
other words, the evolution is 'quasi-unitary'. 



(30) 
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To gain further insight into the dynamics generated by (13"Uj) we note that Q is the 
product of two Bloch spheres. Any motion on Q thus corresponds to a pair of coupled 
trajectories on these spheres. In terms of the usual spherical coordinates ({Oi}, {0i})i=i,2 
a point on Q can be written in the form |V'i)|'02)j where 

= cos§0 i |T>+sin§0 i e 1 *|r). (31) 

The idea of visualising the motions on the spheres is to express the phase-space 
coordinates ({qi}, {Pi}) in terms of the spherical coordinates ({Oi}, {0i})- F° r this, we 
must specify a Hamiltonian so that we can establish the relation between the energy 
eigenstates and the four chosen basis states |Tt)> ITI)? lit)) III) on Q- For example, 
suppose that the Hamiltonian takes the form 

H = -J<ti ® & 2 - B(a z x <g> 1 2 + li <g> o-f), (32) 

which is a Heisenberg-type spin-spin interaction with strength J and an external z-field 
with strength B. The eigenstates of this Hamiltonian are given by the spin-0 singlet 
state and the spin-1 triplet states. By comparing the coefficients of \ipi)\ip2) with (1241) 
one can express ({qi}, {Pi}) in terms of ({Oi}, {4>i})- We invert the resulting relations to 
obtain 

6>i i2 = sin -1 y/p2 +P3- %Jpm, ± cos" 1 y/pi +Pi- y/piPi ^ 
0i,2 = -\ (qi icos" 1 (|(p3 -P2)l^/Pipt)) , 

where (9%, 0i) corresponds to the '+' sign and (62, 02 ) corresponds to the ' — ' sign in the 
right side of fl33l . From fl30l) and fl33|) we deduce the equations of motion in terms of 
the spherical variables, with the results: 

^ = ^ = ° (34) 
01 = 02 = - ^2 - ^3)(sm 2 ±0i + sin 2 ±0 2 ) - \(lo 2 + w 3 ). 

We thus find that for a given initial state the dynamical trajectories are given by a pair 
of latitudinal circles on the respective Bloch spheres. An example of a field plot for one 
of the Bloch spheres is shown in Figure [TJ 

Example 2. We consider a system of three spin-± particles, and impose the 
restriction that the unitary evolution is constrained to lie on a product space 1Z = 
V 1 x V 1 x V 1 that contains the energy eigenstates. In this case, a necessary and 
sufficient condition for the state to lie on 1Z is that the components of the state 
vector simultaneously satisfy two quadratic equations [18]. Expressed in terms of the 
homogeneous coordinates {'0i}i=i,2,...,8 of we nn d that the relevant constraints are 
given by four complex equations: ip\ipj = 1P31P4, ip2ip8 — ^s^s, "01^8 = V^Vv? an d 
1/J31/JQ = ipi^s- Taking the real and the imaginary parts of these equations and expressing 
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the results in terms of ({qi}, {Pi}) we obtain the following eight constraints: 

$ 1 = q 2 - q 5 ~ q& 

$ 2 = qi + q7 ~ <?3 ~ Qa 
$ 3 = qi - q 2 - q 7 

$ = P2P8 - P5P6 

$ 6 = P1 p 7 - pm 

$ 7 = P1P8 - P2P7 
$ 8 = P3P6-P4P5, 

where we have written p$ = 1 — Y^l=iPi- Remarkably, the constraint equations for the 
variables {q{\ and {pt} decouple. We can also read off from fl35l) the condition for the 
unitary motion to lie on 1Z: 

0J\ = 002 + ^7 

u; 3 + lu 6 = uj 4 + to 5 . 

It follows that the eigenvalues of the trace-free part of the Hamiltonian must take the 
form {ei, e 2 , e 3 , e 4 , -e 4 , -e 3 , -e 2 , -ei}. 

If H does not have this property, then the constraints fl35|) become nontrivial. 
Nevertheless, owing to the fact that they decouple it is straightforward to verify that 
the equations of motion are given by qi = fi({pi}) and Pi = for i = 1, . . . , 7, where 
fi({pi}) are elementary functions of the variables {pi}- It follows that the dynamical 
evolution is quasi- unitary in that the amplitudes {pi} remain constant and the relative 
phases {p^} evolve linearly in time. 

Example 3. Let us turn to a different example. We consider again a pair of spin-| 
particles and impose the condition that an initially disentangled quantum state remains 
disentangled under the dynamics. This constraint can be expressed algebraically by 
requiring that the motion is confined to the special quadric Q' that corresponds to 
disentangled spin states. A generic state is still given by ff24l) but it is now no longer 
the case that all the energy eigenstates \Ei) lie on Q'. In other words, the constraint 
surface is no longer given by ipi^A = ^2^3 in the energy basis. 

For the specification of the constraint equation associated with disentangled states, 
we must specify the Hamiltonian. We consider the Heisenberg-type spin-spin interaction 
system given by fl32l) . The corresponding energy eigenstates include the total spin-0 
singlet state and the total spin-1 triplet states. We then change basis by substituting 
the energy eigenstates in (f2lj) with the singlet and triplet spin states such that we can 
express our general state in terms of the four disentangled basis states |||), |f |), ||T)j III) 
of Q'. We find 

m = VWItt) + ^UV2e-' m - vp 3 e~ iqa )\n) 

+ ^(VP2e' iq2 + VP3e' iq3 )\U) + Vpie-^IU). (37) 
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Figure 2. A snapshot of the vector field generated by (0i,<fii) in Example 3. 
The parameters are chosen as £; = |, E 2 = 1, E 3 = —2, E4 = \, and 
O2 = <f>2 = \k. This choice gives the vector field 6\ = cos 4>\ cos 6\ — 3) and 
</>i = I (9 cos 6>i — sin 0i cos 2 Q\ / sin 6*i) . 



Since all the basis states used in ( 1371) lie on Q', now we are able to use the condition 
■01^4 = 4>2'4 ! 3 to constrain the motion on Q': 

VPme~ iqi = |(p 2 e- 2i92 - p 3 e- 2iq:i ). (38) 

Taking the real and imaginary part of (1381) we find that the two constraints required for 
the motion to remain on Q' are 

$1 = Zy/Pm ~ P2 cos(2g 2 - qi) + P3 cos(2g 3 - q{) = 

$2 = p 2 sm(2q 2 - q\) - p 3 sin(2g 3 - g x ) = 0. 1 

Unlike the previous examples, these constraints are not separable, hence we can no 
longer expect to find the quasi-unitary motion seen in the previous cases. Substituting 
the constraints fl39|) into ffT5]) and following the procedures, we find that the resulting 
equations of motion are given by 

pi = 2p 2 p 3 (u 1 - u 3 ) sin(2(g 2 - q 3 )) 
p 2 = -2p 2 p 3 (u)i - 2lj 3 ) sin(2(g 2 - q 3 )) 
p 3 = 2p 2 p 3 (uj 1 - 2^2) sin(2(g 2 - q 3 )) 

q 1 = 2p 3 (l- 2p 1 -p 2 -p 3 )(uj 2 - uj 3 ) cos(2g 3 - qi)/^/pm 

+(wi - 2oo 2 )(2 Pl + p 2 + p 3 ) + 2u 2 
q 2 = 2p 1 p 3 (u 3 - u 2 ) cos(2g 3 - qi)/^/p~m + (wi - 2u 2 )(2p 1 + p 2 ) -p 3 (^i - 2^ 3 ) + 2u 2 
q 3 = 2p 1 p 3 (u 3 - u 2 ) cos(2g 3 - q{)/ y/pip* 

- 2^2) (2pi -p 2 cos(2(q 2 - q 3 ))) + p 3 (uJi - 2u 3 ) + 2^ 3 . 

It is no longer the case that pi = 0. Therefore, we have six coupled nonlinear differential 
equations describing the motion. We can visualise the motion on the product of two 



Constrained quantum motion 



11 



Bloch spheres by converting f)40p into spherical coordinates: 

9\ = sin(0i - 02 ) sin^ 2 [(^i - ^2) cos^i + uj 2 - iv 3 ] 
9 2 = sin(0i - 02 ) sin 9 X [(^2 - ^1) cos 9 2 - tu 2 + lu 3 ] 

01 = \[-uji + (co> 2 - ^) cos ^ 2 + (f^i - ^2 - 2cj 3 ) cos #1 + [cos(0i - 2 )/(sin6 l i sinfto)Jx 

x (2(o; 3 — uo 2 ) sin 2 #1 cos# 2 + (^1 — ^ 2 )(cos 2 9\ — cos 2 9 2 ))] 

02 = + (cu 2 - ^) cos6q + (|o;i - uj 2 - 2cj 3 ) cos# 2 + [cos(0i - 2 )/(sin0i sin6> 2 )] 

x (2(0*3 — o> 2 ) cos^i sin 2 9 2 — (&i — a;2)(cos 2 9i — cos 2 ^2))] • 

An example of the resulting field plot arising from these equations is shown in Figure [21 
indicating the nontrivial nature of the dynamics. 

In summary, we see that constrained quantum motions of the Dirac type can be 
treated straightforwardly by means of the prescription described above. The resulting 
dynamical equations are in general quite intricate, and typically require numerical 
analysis. 

The authors thank D. Holm and M. Parry for stimulating discussions. This work was carried out 
in part while DCB and LPH were visitors at the Perimeter Institute for Theoretical Physics, Waterloo, 
Ontario. 
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